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Abstract 

We calculate the leading order QCD fragmentation functions for the production of 
P-wave charmed beauty mesons. Long-distance effects are factored into two nonpertur- 
bative parameters: the derivative of the radial wavefunction at the origin and a second 
parameter related to the probability for a (be) heavy quark pair that is produced in 
a color-octet 5"-wave state to form a color-singlet P-wave bound state. The four 2P 
states and those 3P states which lie below the BD flavor threshold eventually all decay 
into the IS ground state B c through hadronic cascades or by emitting photons. The 
total fragmentation probabilities for production of the 15" ground state B c from the 
cascades of the IP and 3P states are about 1.7 x 10~ 4 and 2.3 x 10 -4 respectively. 
Thus the direct production of the P-wave states via fragmentation may account for a 
significant fraction of the inclusive production rate of the B c at large transverse mo- 
mentum in high energy colliders. Our analytic results for the P-wave fragmentation 
functions disagree with those obtained earlier in the literature. Our results also cover 
the P-wave heavy quarkonium case in the equal mass limit. 



1 Introduction 



Much progress has been made in the past several years to improve our theoretical under- 
standing of the physics of a heavy hadron containing a single heavy quark (charm or bottom 
quark) . This is mainly due to the discover of the powerful heavy quark spin- flavor symmetry 
[0] in Quantum Chromodynamics (QCD) and the development of the Heavy Quark Effective 
Theory 0. Heavy hadrons containing both the heavy bottom and charm quarks are also 
of interest. Of particular importance is the (be) meson family. Since the top quark is very 
heavy and will decay rapidly into W + b before hadronization, the (be) system might be the 
only physical system containing two heavy flavors of different masses. Techniques of Heavy 
Quark Effective Theory can not be applied here because the charm quark is not sufficient 
light. However, since this family is intermediate between the cc charmonium and the bb bot- 
tomonium systems, many properties (mass spectrum, transition and decay rates, lifetimes, 
etc.) of this physical system can be studied by using potential models and QCD sum rules 

ill- 

Recently, it has been pointed out || that the dominant production mechanism for 
these yet unobserved (be) bound states with large transverse momentum at a high energy 
process is fragmentation, in which an energetic b antiquark is first produced at large trans- 
verse momentum by a short-distance hard process and subsequently fragments into various 
(be) bound states. Furthermore, the process independent fragmentation functions describing 
these phenomena are shown to be calculable at the heavy quark mass scale using perturba- 
tive QCD. In particular, the fragmentation functions for the b antiquark to split into the two 
S'-wave states ^(rz 1 ^) and B*(n 3 Si) were calculated |5j, |6| to leading order in both a s and 
v, where v is the typical relative velocity of the charm quark inside the meson. Including the 
IS and 2S states that lie below the BD flavor threshold (M threshold = M D + M B « 7.1 GeV), 
a lower bound on the inclusive branching fraction for the production of the B c ground state 
from Z° decay has been estimated to be about 2.3 x 1CT 4 ||, which might be too small for 
this particle to be observed at LEP with the present luminosity. Production of the B c meson 
at the hadron colliders via the direct b antiquark and induced gluon fragmentation has also 
been studied in Ref.[[7[. About 20,000 B c with transverse momentum > 10 GeV are 
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expected to be produced at the Tevatron with an integrated luminosity of 25 pb . A clean 
signature of B c would be the observation of three charged leptons coming from the same 
secondary vertex due to the decay — > J/ip + l'vi< followed by J/ip — > U- The combined 
branching fraction for these decays is expected to be about 0.2% ||. One expects about 40 
of these distinctive events at the Tevatron. Thus, unless LEP can increase its luminosity 
by an order of magnitude or so in the near future, the best place to look for the charmed 
beauty meson will be at the Tevatron. 

The mass spectrum of the (be) mesons has been obtained by Eichten and Quigg || 
using potential model calculations. Since, unlike the quarkonium system, the b and c quarks 
do not annihilate into gluons, the excited states of the (be) mesons can not decay directly 
into light hadrons. States that lie above the BD threshold will strongly decay into a pair of 
B and D mesons. States that lie below the BD threshold will decay eventually to the ground 
state B c either by hadronic cascades or by emitting photons. In order to get more reliable 
theoretical prediction for the inclusive production rate of the B c meson at large transverse 
momentum, one must include all higher orbitally excited states of the (be) mesons that lie 
below the BD threshold. Below the BD flavor threshold, there are two sets of S-w&ve states 
(the IS" and 2S states), as many as two complete sets of P-wave states (the 2P and perhaps 
some or all 3P states), and one set of P-wave states (the 3D states). 

In this paper, we report the results of the calculation of the fragmentation functions 
for b antiquark splittings into the four P-wave (be) mesons - the X P\ and 3 Pj (J = 0, 1,2) 
states. In Section 2, we write down new factorization formulas for the P-wave fragmentation 
functions that are valid to all orders in strong coupling constant a s and to leading order in 
v, where v is the typical relative velocity of the heavy quarks inside the (be) meson. Long- 
distance effects are factored into two nonperturbative parameters: the derivative of the radial 
wavefunction at the origin, and a second parameter related to the probability for the (bc) pair 
that is produced in a color-octet S"-wave state to form a color-singlet P-wave bound state. In 
Section 3, we calculate the perturbative coefficients of the fragmentation functions to leading 
order in a s . Infrared logarithmically-divergent terms that appeared at the next-to- leading 
order in a s are isolated and absorbed into the second nonperturbative parameter. Mixing 
effects between the l Pi and 3 Pi states in the fragmentation functions are also calculated. In 
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Section 4, we take the heavy quark limit of mb/m c — > oo, and show that our fragmentation 
functions exhibit the heavy quark spin symmetry while earlier results obtained by Chen || do 
not satisfy this important property. In Section 5, we discuss a crossing relation between the 
fragmentation function D^h(z) for parton i splittings into hadron H and the distribution 
function fi/niz) of finding parton i inside hadron H. In Section 6, we consider the equal 
mass limit corresponding to the case of P-wave heavy quarkonium. Numerical results and 
conclusions are made in Section 7. 



2 New factorization formulas for P-wave fragmenta- 
tion functions 

Recently Bodwin, Braaten, and Lepage JTO] have reformulated the perturbative calculations 
of the annihilation rates of heavy quarkonium using the framework of non-relativistic QCD 
(NRQCD) jl0[ [lip. This work points out that the usual factorization assumption that all 
long-distance effects can be absorbed into the radial wavefunctions at the origin and their 
derivatives is not correct when higher order QCD corrections and/or relativistic corrections 
are taken into account. Infrared divergences that arise in these higher order calculations spoil 
factorization. To obtain sensible perturbative results that are free from infrared divergences, 
one must use new factorization formulas that involve a double expansion in the strong cou- 
pling constant a s and in the typical velocity v of the heavy quark inside the quarkonium 
[ |10|j . The easiest way to organize the expansion in v in these theoretical calculations is to 
use NRQCD 



As in the quarkonium case |12[], the heavy quarks inside the (be) meson are mov- 
ing nonrelativistically and separated at a typical distance of order l/(mv), where v -C 1 is 
the typical relative velocity of the heavy quarks and m = m^mc/ (m.& + m c ) is the reduced 
mass. The calculation of the perturbative fragmentation function is based on separating 
short-distance perturbative effects involving the scale of order 1/m from long-distance non- 
perturbative effects involving scales of order l/(mv ) or larger. To leading order in v , two dis- 
tinct mechanisms that contribute to the fragmentation function have been identified and are 
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referred as the color-singlet mechanism and color-octet mechanism. The color-singlet mech- 
anism is the production of a (be) pair in a color-singlet 1 P\ or 3 Pj state with separation of 
order 1/m in the (6c) meson rest frame. The subsequent formation of the P-wave (be) bound 
state is a long-distance nonperturbative process with a probability of order v 2 (mv/m) 3 ~ v 5 , 

where the first factor v 2 arises from the derivative of the wavefunction of the P-wave state 
and the second factor (mv/m) 3 is a volume factor |I3J. The color-octet mechanism is the 

production of the (be) pair in a color-octet x Sq or 3 Si state with a separation of order 1/m. 
The subsequent formation of the P-wave (be) bound state can proceed either through the 
dominant \bc) component or through the small \bcg) component of the wavefunction. In the 
first case, the (be) pair must radiate a soft gluon to make a transition to the color-singlet 
P-wave bound state. In the second case, a soft gluon must combine with the (be) pair to 
form the color-singlet P-wave bound state. In both cases, the probability is of order v 5 , with 
a factor of v 3 coming from the volume factor and an addition factor of v 2 coming either from 
the probability of radiating a soft gluon or from the small probability of the \bcg) component 
of the wavefunction pi 113]. Since the color-singlet mechanism and color-octet mechanism 



contribute to the fragmentation function at the same order in v , they must both be included 
for a consistent calculation. 

Separation of the short- distance physics from the long-distance effects discussed in 
the previous paragraph requires the introduction of an arbitrary factorization scale A in 
the range mv C ACm. The fragmentation functions for the P-wave (be) mesons satisfy 
factorization formulas that involve this arbitrary scale. To leading order in v 2 and to all 
order in a s , the factorization formulas for the fragmentation functions of b —>■ bc(n l Pi) and 
b — > bc(n 3 Pj) at a scale /io near the heavy quark mass scale consist of two terms: 



D (■ // ) H Kbc){n) (1) #8(6c)( A ) D ( 8 ) , 



D^ Hn s Pj) (z,» ) = ^^DM- bc{3pj) ^ (2) 



where Hi and H' 8 (A) are nonperturbative long-distance factors associated with the color- 
singlet and color-octet mechanism respectively, m is the reduced mass, and n — 2, 3, • • •, labels 
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the principal quantum number of the P-wave states. The short- distance factors D^(z,A) 
and D^(z) can be calculated perturbatively as power series in a s . They are proportional to 
the fragmentation functions for a b antiquark to split into a (be) pair with vanishing relative 
momentum and definite color-spin-orbital quantum numbers: color-singlet 1 P\ or 3 Pj state 
for and color-octet 1 So or 3 Si state for D^ 8 \ Note that in these factorization formulas, 
the only dependence on A is in and H' 8 . As in the quarkonium case, the nonperturbative 
parameters Hi and H' 8 can be rigorously defined as matrix elements of 4-quark operators in 



nonrelativistic QCD divided by appropriate powers of the heavy quark masses ||10fl . Their 
dependence on A is governed by a renormalization group equation whose coefficient can be 
calculated using nonrelativistic QCD [TO]. We adopt the same definitions for these matrix 
elements as in the charmonium case discussed in Ref.|lO|, except that explicit factors of m c 



are replaced by 2m. To order a s , H 1 Q )C - ) (n) is scale invariant and is related to the derivative 
of the non-relativistic radial wavefunction at the origin for the P-wave (be) bound states: 

Jf%1? ( i+o( " 2) ) ; (3) 



while HLj^ satisfies |T(], [14 



A ^ H 'sr bc M) = ^-a s (A)Hu ic) , (4) 



with the solution 



Kile) (A) = ^ ( - 6c) (Ao) + ^log ) H l{lc) , (5) 



where /3q = 9/2 is the first coefficient in the beta function for QCD with 3 flavors of light 
quarks. If the factorization scale A is chosen to be much less than the reduced mass m, the 
above equation can be used to sum up large logarithms of m/A. To avoid large logarithms 
in the perturbative coefficients D^ l > in Eqns.(|I]) and (Q), one can choose A on the order 
of m. Since |P^p(0)| 2 is of order v 5 |i~0| , from Eqns.(|3]) and we see that the color- 



singlet and color-octet contributions in the fragmentation functions are both of order v 5 . 
Following Ref . p[ , the initial scale fiQ for the P-wave fragmentation functions in Eqns. (ffl) and 



(0) can be chosen to be (m& + 2m c ) - the minimal virtuality of the fragmenting b antiquark. 
Fragmentation functions at a higher scale can be obtained by solving the Altarelli-Parisi 
evolution equation with Eqns.([3]) and (@) as the boundary conditions. These boundary 
conditions will be determined in the next Section. 



3 Fragmentation functions for b — > P-wave (6c) mesons 

We now embark on the calculation of the coefficient D^(z,A) in the color-singlet con- 
tribution to the fragmentation function. Let us first briefly review the method that was 
introduced in Refs.[[5], [15[ of calculating the heavy quark fragmentation functions using per- 
turbative QCD. We refer to Ref.|| for more details. Let M. denote the amplitude for a high 
energy source (symbolically denoted by T) to create a b*(q) — > H(p) + c{p') with total 4- 
momentum q = p+p' as illustrated in Fig.l. Here H denotes a (be) bound state. The leading 
order diagram of Fig.l involves creating a (cc) pair from the vacuum. Let .Mo denote the 
amplitude for the same source to create an on-shell b antiquark with the same 3-momentum 
q. Then the fragmentation function for b — > H is given by || 

n ( s 1 f , J {m b + m c f m 2 c \ \M\ 2 

D- b ^ H {z) = ds6[s hm — — , 6 

167?^ J \ Z 1 — Z j q /m b , c ^oo \Mo\ 

where s = q 2 is the virtuality of the fragmenting b antiquark. In a frame where the virtual 
b antiquark has 4-momentum q = (q , 0, 0, (73), the longitudinal momentum fraction of H is 
z = {Po + P3)/(?o + ^3) an d its transverse momentum is p± = (pi,^)- The relation between 
the transverse momentum and the virtuality is given by 

^ = z(l-^-_-_j. (7) 

The tree-level matrix element squared |.Mo| 2 is simply given by 3Tr(rf^) as qo/rrit, — > 00. 
The matrix element M. for producing a bound state H(p) consisting of collinear b and c 
quarks with relative momentum k (Fig.l) can be expressed as a Dirac trace involving the 
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spinor factor v(rp + k)u(rp — k). To project out the (be) bound state in the non-relativistic 
approximation, this spinor factor is replaced by the following projection operator 



i(rp + k)u(rp — k) — > y/rrib + 



' r ft - m b \ (r p 1 — f + m c 



2rrih 



2m r 



(8) 



with O = 75 or -£*(p, S z ) for the spin singlet or spin triplet state respectively, where e(p, S z ) 
is the polarization vector associated with the spin triplet state. Besides this spinor factor, 
the gluon propagator in Fig.l also depends on the relative momentum k. In the axial gauge 
associated with the 4-vector n^, the gluon propagator is given by 



1 



- g» v + ; — + 



n ■ I 



(n ■ If 



(9) 



where / = q — fp — k. The axial gauge vector can be chosen to be = (1, —p/\p\) so that 
Fig.l is the dominant diagram and hence factorization is manifest. By using the equations 
of motion for the b and c quarks, one can show that terms that are proportional to l v in 
the gluon propagator do not contribute to the matrix element M.. For the P-wave states, 
one needs to expand the above projection operator and the gluon propagator appearing in 
the matrix element M. to first order in the relative momentum k. After some manipulations 
by using the standard covariant methods [16], the matrix element for the n^Pi state can be 
written in the following form: 



M(n x P x ) =5, 



1 g 2 R' nP (Q) 



1 



lJ 3r 2 f y^M (s-f 2 M 2 ) 3 



e* a (p,L z )TV a l5 v(p') 



(10) 



where % and j are the fundamental color indices, g is the strong coupling constant, and 
e Q (p, L z ) is the polarization vector of the 1 P\ state. The vertex factor V a is given by 



ya = Sfq a (d+fM)(p'-2M) 



- (s-f 2 M 2 

(s - f 2 M 2 ) 2 
n ■ (q — fp) 



2(1 -2r){4+ fM) 7 a - 



4r 



n- (q — rp) 



q a (p'+M)i 



— (1 - 2r)M) 7 a - 2rf — — —n a (p' + M) 

M n ■ (q — rp) 



i- (ii) 
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In Eqn.(|TT|), we have denned M = + m c , r = m c /M, and r = 1 — r = rrib/M. Similarly, 
the matrix elements M. (n 3 Pj) for the three spin triplet states can be written in the following 
compact form: 



M(n 3 Pj) = 5, 



1 g 2 R' nP (0) 



1 



v 3r 2 f jAnM (s-f 2 M 2 ) 3 



^ a p{J,J z )TV^v{p'), 



(12) 



with 



*«fi(J,Jz) = J2(l,L z ;l,S z \J,J z )e* a (p,L z )e* p (p,S z 

^Sap-W) forJ = 0, 

-^iftctfcvP* 671 * (P) J z) for J = 1 , 
tapijP, Jz) for J = 2, 



(13) 



(14) 



where e a (p, L z ) and e^p, S z ) are the polarization vectors in the orbital and spin space of the 
3 Pj states respectively; and e a (p, J z ) and e a p(p, J z ) are the helicity wave functions of the total 
spin J = 1 and J = 2 states respectively. Our conventions are g^ u = diag(l, — 1, — 1, — 1) 
and e 0123 = +1. The vertex tensor V a @ in Eqn.(|12|) is given by 



V 



af3 



lfMq a {q' + fM)~/ /3 



+ (s-f 2 M 2 ) 



^U+fM) (2^((1 - 2r) i- 2M) + [ 7 *,/] ^) 



+ 



+ 



4r 



n ■ (q — rp) 

(s - f 2 M 2 ) 2 
n ■ (q — fp) 

1 



—q a (p'+M) 1 p fi 



-L (^((1 - 2r) jrf + M) - (1 - 2r)M)[ 7 a , 7 ^] 



n ■ (q — rp) 



i 



(15) 



The general procedure for extracting the P-wave fragmentation functions from the 
above matrix elements is the same as in the case of the .S-wave || except the intermediate 
steps are much more involved. We will skip over these tedious steps and present only the 
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final results. From Eqn.fQ), we can define generalized fragmentation functions D^(z, s) that 
depend on both the longitudinal momentum fraction z and the virtuality s according to 



M- r-'.\/- 



z 1 — z ) 



(16) 



For the x Pi state, we have 



with 



81 (l-f^) 4 ^ (s-f 2 M 2 ) 5 -' ; 



(17) 



/o = 64r 2 r 3 (l — rz 



/i = 8tt(1 — rz) 



h = -iX-fz) 



3 - 2r - 2r - 2f (2 + 4r - r 2 )z + f 2 (l - 2r)z 



2(1 - 2r + 4r 2 ) - (3 - 42r + 64r 2 - 16r 3 )z 



(18) 
(19) 



2rf(23 - 14r - Ar 2 )z 2 + f\l + 12r)(l - 2r)z d 



(20) 



/3 = (!"*) 



1 - 2(1 - 2r)z + (3 - 2r + 2r 2 )z 2 - 2f(2 + r - 2r 2 )z 3 + f 3 (2 + r 2 )z 4 



(21) 



Integrating Eqn.(|l7|) over s as in Eqn.([l6|), we obtain the fragmentation function D^(z) for 
the color-singlet contribution in the 1 P\ state 



D^- (z) 



16a 2 (2m c ) rr 3 ;z(l — z)' 



{6-6(4r 2 -8r + 5); 



243 (1 - fzf 
+ (32r 4 - 96r 3 + 250r 2 - 210r + 69)z 2 
+ 8f(4r 4 + 12r 3 - 48r 2 + 37r - 12)z 3 

+ 2f 2 (16r 4 + 161r 2 - 114r + 42 > 4 - 6f 3 (4r 3 + 28r 2 - 15r + 7)z 5 

+ f 4 (46r 2 - 14r + 9)z 6 } . (22) 
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Similarly, for the 3 Pj states, we have 



D 



(i) 



(z, s) 



32a 2 (2m c ) rf 



3 3 



E 



2n 



243 (1 - fzY ^ {s - f 2 M 2 f~ n ' 



with 



/ <°> = 64r 2 f 3 (l - fzf , 



f[ ^ = 8rr(l — rz) 1 



1 - 18r + 14r 2 - 2f(l - 2r + 7r 2 )^ + f 2 (l + 2r)z 2 



/f = -(l-f*) s 



2(1 - 4r)(l + 6r - 4r 2 ) - (5 + 14r - 8r 2 + 80r 3 - 64r 4 )z 



+2f(2 + 9r + 18r 2 - 28r 3 - 16r A )z 2 - f 2 (l + 6r + 16r 2 - 32r 3 )/ 



(o) 



(1-z) 



1 — 4r — (1 — 4r)(l - 2r)z - rr(3 - 4r)^ 2 



/ (1) = 192r 2 f 3 (l - fzf , 



= 24rf(l-f^) ; 



2(1 - r - r 2 ) - f(3 + lOr - 2r 2 )^ + f 2 z 2 



ff ] = -6(1 -fzf 



2(1 + 2r) - (5-2r + 6r 2 )z 



+2f(2 - 3r - 4r 2 )^ 2 - f 2 (l - 2r + 2r 2 )^ 3 



/ 3 (1) = 6(1 -z) 



1 - 2(1 - 2r)z + (1 - 4r)(l - 2r)z 2 + 2rf(l - 2r)^ 3 + r 2 f 2 z 4 



and 



/ {2) = 320r 2 f 3 (l-fz) 4 , 



/f) = 8rf 2 (l-f^ 



2(4 + 13r) - (1 + 70r - 26r 2 )z - f(7 + 8r)^ 2 



/ 2 (2) = -^(l-fzf 



4(1 + 4r) - (7+ 12r-32r 2 )^ 
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+2(1 + 13r - 26r 2 + 8r 3 )/ + (1 - 30r - 5r 2 + 4r 3 )z 



3\ J 



(34) 



/f = 4f 2 (l-z) 



2 - 4(1 - 2r)z + (5 - 8r + 12r 2 )/ 



2(1 - 2r)(3 + 2r 2 )/ + (3 - 12r + 12r 2 + 2r 4 )/ 



(35) 



Integrating Eqn.(|23"D over s as in Eqn.(|T6|), we obtain 



D 



(i) 

fe^fec(3P ) 



16a 2 (2m c ) rr 3 z(l — z) 2 ( , ; 9 



729 (1 - f Z y 
+ (832r 4 - 1456r 3 + 1058r 2 - 362r + 63)/ 

- 8f(100r 4 -184r 3 + 118r 2 -22r + 9)/ 
+ 2f 2 (416r 4 - 776r 3 + 369r 2 + 42r + 24)/ 

- 2f 3 (240r 4 - 516r 3 + 232r 2 + 59r + 9)/ 
+ f 4 (96r 4 - 240r 3 + 134r 2 + 34r + 3)/} , 



(36) 



and 



32a 2 (2m c ) rf 3 ^(l - z) 
243 (1 - fzf 



{6 + 6(4r-5)z 



+ (16r 4 + 64r 2 - 98r + 63)/ + 8f(2r 4 + 2r 3 - 13r 2 + llr - 9)/ 

+ 2f 2 (8r 4 - 16r 3 + 47r 2 - 18r + 24)/ + 2f 3 (8r 3 - 24r 2 + r - 9)z 5 

+ f 4 (12r 2 + 2r + 3)/} , (37) 



-Dfe^ bc (3p 2 )(^) 



64a 2 (2m c ) rr 5 z(l — zf 



{12 + 12(2r-5)z 



729 (1 - fz) 
+ (92r 2 - 76r + 135)/ + 4(10r 3 - 54r 2 + 31r - 45)/ 
+ 2(46r 4 - 16r 3 + 123r 2 - 78r + 75)/ 
- 4f(6r 4 + 9r 3 + 40r 2 - 13r + 18)/ 
+ f 2 (12r 4 - 12r 3 + 55r 2 - lOr + 15)/} . 



(38) 
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As in Ref. H], the scale of the strong coupling constant is chosen to be 2m c - the minimal 
energy of the exchange gluon. The results of these fragmentation functions for the P-wave 
states disagrees with those given in Ref.[(|. 



Integrating Eqns.(|2~2|) and (|36"D-(|38"D over z, we obtain the fragmentation probabilities 
P^ = (Hi/m) Jq 1 dz D^(z) of the color-singlet contribution for b fragments into bc(n 2S+1 Pj): 



P-, 



(i) 



16a 2 (2m c ) H 1( j, c) (n) 



243 



ni 



!^f^(10r 4 + 50r 3 + 91r 2 + 20r + 7) 



35r 



r (24r 6 - 256r 5 - 2083r 4 - 9538r 3 - 5758r 2 - 907r - 172) 



(39) 



P (i) 

b^bc(n 3 P ) 



16a 2 (2m c ) H m (n) 



729 



r logr 



rn 



3— ^(32r 6 + 48r 5 - 398r 4 + 194r 3 + 337r 2 + 32r + 1) 



: (9944r 6 - 17384r 5 - 34289r 4 + 56116r 3 + 11036r 2 + 347r + 60) 



35r : 



(40) 



6->6c(n 3 Pi) 



32a 2 (2m c ) H 1(ic) (n) 



243 



rn 



!^^(20r 4 + 46r 3 + 31r 2 + 8r + 1) 



35f 5 



;i2r 6 - 268r 5 - 4461r 4 - 4006r 3 - 2326r 2 - 23r - 58) 



(41) 



and 



Pi 



(i) 



64a 2 (2m c ) H m {n 



b^bc(n 3 P 2 ) 



729 



rn 



3 rlogr ^ 6 + 3g ^ 5 + ^ 4 _ 20 ^ 3 + 140 ^ 2 + Ar + n ^ 



+ — -(1468r 6 + 5570r 5 - 7345r 4 + 12080r 3 + 6130r 2 + 712r + 285) 
35r 2 



(42) 
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For a given principal quantum number n, the 1 P± state and the 3 Pi state, constructed 
in the LS coupling scheme, are mixed in general to give rise to the physical states 1 + and 
1 + ' defined as 



|1+') = cos^Pi) - sin0| 3 Pi) , 
|1 + ) = sin^Pi) + cos0| 3 Pi) . 



(43) 
(44) 



Thus, in general, we have 



D 



(i) 
i— +/ 

(i) 



D lZ H i + ')( z ) = cos 2 6D^- bc{lpi) (z) + sin 2 6D^ (45) 



'5-fc(i+)W sin^P-^ (1Pi) ( 2 ) + cos^P^ c(3Pi) (z) + sin^cos^4W- (46) 

Similar equations hold for the corresponding generalized fragmentation functions D^'(z, s) 
as well. With the matrix elements of M.(n l P\) and Ai(n 3 Pi) given by Eqns. (|T0|) - (|11|) and 
(|T2p - ( ITH ) respectively, we can calculate their interferences and obtain the mixing term: 

64 v / 2a 2 (2m c 



(i) 



tyiy 3 



81 



2 ff x M 6 ~ 2n 



(l-rz) 3 ^ (s-f 2 M 2 ) 



4— n 



with 



mix 




]rr(l — rz) 3 



2 — 3r — r(l + r)<z 



J I 



[l-fz] 



f2 



'l-z) 



2(1 + 2r 2 ) - (5 - 24r + 36r 2 - 8r 3 )z 



+2f(2 - 15r + 10r 2 - 2r 3 )^ 2 - f 2 (l -Br- Ar 2 )z 3 



1 - (1 - 3r)z + r(2 + r)z 2 + frV 



(47) 



(48) 



(49) 



(50) 



Integrating Eqn. fl47D over s as in Eqn. fllq) , we obtain 



32v / 2a 2 (2m c ) rf 3 z{l-zf 



{6-6(2r 2 -4r + 3)^ 



243 (1 - fzf 

+ (24r 3 + 52r 2 - 52r + 21)z 2 + 2f(14r 3 - 6r 2 + 15r - Q)z 3 



r 2 



(2r 2 + 8r - 3)^ 4 } . 



(51) 
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The fragmentation probability of the mixing term is given by 



mix 



(1) 



r i(bc){ n ) f 1 r n (l) / s 

m Jo 

32v / 2«^(2m c ) H 1( T. c) (n) \ rlogr, , ,, 9 

— ^ c ± 1(bc)V ' 3 — |- 18r 4 + 28r 3 - llr 2 - 2r + 1) 

243 m f 3 v y 



+ ^r(24r 5 + 587r 4 - 18r 3 - 88r 2 - 3r + 8) 



(52) 



We next turn to the calculation of the color-octet coefficient D^(z) to leading order 
in a s . Notice that to leading order in a s , all the color-singlet contributions obtained above 
are free from infrared divergences. All the long-distance effects can be factored and absorbed 
into the derivative of the wavefunction, or equivalently Hi to leading order in v 2 . The color- 
singlet coefficients D^(z) do not depend on the factorization scale A to leading order in a s . 
Beyond leading order, infrared singularities appear in the color-singlet contribution arise from 
the radiation of a soft gluon in the process b* — > (be) eg. Imagine we are actually doing such a 
next-to-leading order calculation of D^(z) with an energy cutoff A (mv C A<m) imposed 
on the radiated soft gluon. This cutoff will allow us to isolate the infrared logarithmic 
divergent terms in this higher order calculation and absorb them into the nonperturbative 
parameter H' 8 associated with the color-octet mechanism. However, there is a short cut to 
achieve the same goal by following the same method as in Ref.||12||. One can calculate the 
fragmentation functions for the processes b* — > &c( 1 5 , ,8) and b* — > bc( 3 Si,8), where the (be) 
pair is in the appropriate color-octet S"-wave state. The calculations of these fragmentation 
functions are identical to the color-singlet S'-wave case ||, except one has to replace the 
S- wave wavefunction Rs(0) by a fictitious "color-octet S'-wave wavefunction at the origin" 
-Rs(O), and the color factor C\ by CV/6 where Cj? = 4/3 for color SU(3). Rs(0) is related 
to the nonperturbative parameter HLtJA) by 



Using this trick, we can easily extract the coefficients of the color-octet contributions D^ 8 \z). 



#8(5e)(A) 



2 |i? 8 (0)| 2 
3vr (2m) 2 



(53) 
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Analogous to Eqn.(|16"D of the color-singlet case, we define D^ 8 '(z,s) as the following: 

D®{z) = JdseL-^-- r -^jD^\z, s) , (54) 



where 



£)(8)( 



z, s 



a 2 (2m c ) rr 3 



2 o„M 6 " 2n 



E 



1 (1 - rz) 2 ^ (s - f 2 M 2 ) A - n 



(55) 



For the 1 5'o state, we have 



S'o(iSo) = -12rr(l - rz) 



0i( l s o ) = ~(l-rz) 



2(1 + 2r) - (1 + 12r - 4r 2 )z - f(l + 2r)z 2 



1 + 2rz + (2 + r 2 )z 2 



(56) 



(57) 



(5f 



Similarly, for the 3 Si state, we have 

— 4rf(l — fz) 2 



#0(3 so 



9i(*Si) = -(1-rz) 



2(1 - 2r) - (3 - 4r + 4r 2 )z + f(l - 2r)z 2 



02(BSi) = (l-z)(l+rz)' 



(59) 



(60) 



(61) 



Integrating Eqn. Q55|) over s as in Eqn.(|5"4T), we obtain the color-octet contributions to the 
fragmentation functions 



a 2 (2m c ) rr 3 z(l — z) £ 
162 (1 - fz) 6 



6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )z 2 



- 2f(6 - 19r + 18r 2 )z 3 + 3f 2 (l - 2r + 2r 2 )z 



(62) 



and 



a 2 (2m c ) rr 3 z(l — z)' 
162 (1 - fz) 6 



2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )z 2 



- 2f(4 - r + 2r 2 )z 3 + f 2 (3 - 2r + 2r 2 )z 



„2\„4 



(63) 
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Integrating Eqns.(|62"D and ( |63"D over z, we obtain the fragmentation probabilities P^ 
(ifg/m) Jq dz D^\z) for the color-octet contributions: 



3(8) 

b->bc(}So) 



a 2 s {2m c )H' m (k) 



54 



m 



r logr 



1 + 8r + r 2 - 6r 3 + 2r 4 ) 



+ 



15r 



+ 13r + 228r 2 - 212r 3 + 53r 4 



(64) 



and 



d(8) 
fe^6c( 3 Si) 



a 2 (2m c )^ c) (A) 



162 



m 



r logr 



(7 - 4r + 3r 2 + 10r 3 + 2r 4 



+ -—7(24 + 109r - 126r 2 + 174r 3 + 89r 4 
\hr z 



(65) 



To avoid large logarithms of to/A appearing in the higher order calculation of the color- 
singlet coefficients D^'(z), one should choose A ~ m in the matrix element HL^JA). 

Combining both contributions from the color- singlet and color-octet mechanisms, and 
including the mixing effects in the 1 P\ and 3 Pi states, the total fragmentation probabilities 
are given by 



p-b- 


^bc(n 3 P ) 


p(l) , p(8) 
• r &-*6c(ri 3 .Fb) ^ft-^c^Si) ' 






(66) 


p-b- 


»&c(nl+') 


= cos 2 ^ c(nlPi) + sin 2 ^ c(n3Pl) 


— sin 6 1 cos ^-Pmix 


4- 3 P-®- 


(67) 


p-b- 


->bc(nl+) 


= sm 2 ^5 bc(nlPi) + cos 2 ^ c(n3Pl) 


+ sin 6 1 cos 9Pjri r 


4- 3 P-®- 


(68) 


p-b- 


+6c(n 3 P 2 ) 


p(l) +c:p(8) 
— ^b^bc(n 3 P 2 ) " r ,J - r b^5c( 3 S , i) • 






(69) 



As in the S- wave case 0, under the Altarelli-Parisi evolution to a higher scale, the shapes 
of these P-wave fragmentation functions are softened while the fragmentation probabilities 
remain unchanged. 

Before closing this Section, we note that one can also define fragmentation functions 
that depend on both the longitudinal momentum fraction z and the transverse momentum 
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\p±\- Introducing the dimensionless variable t = \p±\/M and using Eqn.(^) to trade the 
variable s with t, we can define the generalized fragmentation functions D(z,t) and D(t) 
according to 



/ dt D[t) = / dt dz D(z,t) 
Jo Jo Jo 



(Is I dz ois 



r 2 M 2 
1-z 



\D(z,s) 



(70) 
(71) 



Therefore, 



D{z,t) 



2MH 
z(l-z 



■D(z, s) , with s 



M 



1 + t 2 



+ 



r 2 + t 2 



(72) 



The above relation holds for the color-singlet and the color-octet contributions. These func- 
tions D(t) and D(z, t) are useful for studying the transverse motion of the meson with respect 
to the jet axis of the fragmenting heavy quark. The expressions of D(t) for the color-octet 
S- wave contributions can be obtained from Ref.[T^]. Analytic results of D(t) for the P-wave 
states can also be derived but will not be given here. 



4 Heavy quark symmetry 

In this Section, we temporarily leave the heavy-heavy (be) system and turn our attention to 
the heavy- light B- or D-meson system. In the limit of mq/ Kqcd — > oo, both the heavy quark 
spin Sq and the total spin J of a heavy hadron containing a single heavy quark Q become 
good quantum numbers. This implies that in the spectroscopy of the hadron containing a 
single heavy quark Q, the angular momentum of the light degrees of freedom J\ = J — Sq 
is also a good quantum number. We refer collectively to all the degrees of freedom in the 
heavy-light hadron other than the heavy quark as the light degrees of freedom. For the 

heavy-light (Qq) meson, Ji = S q + L where S q is the spin of the light antiquark q and L is 
the orbital angular momentum. Thus the hadronic states can be labeled simultaneously by 
the eigenvalues j and ji of the total spin J and the angular momentum of the light degrees 



of freedom J\ respectively. In general [18], the spectrum of hadrons containing a single 
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heavy quark Q has, for each ji, a degenerate doublet with total spins j + — ji + 1/2 and 
j_ = ji — 1/2. (For the case of jz = 0, the total spin must be 1/2.) For P-wave heavy-light 
mesons, the orbital angular momentum L = 1 and ji can then be either 1/2 or 3/2. Thus 
= (0, 1) and (1,2) for j t = 1/2 and 3/2 respectively. As a result, we expect to have 
two distinct doublets ( 3 Po, 1 + ) and (1 + , 3 P2) in the limit of mf,/m c — > 00, i.e. r — > 0. In this 
limit, the mixing coefficients in Eqns.([43"D-(f44D can be determined by the Clebsch- Gordon 
coefficients in the tensor product of a spin 1/2 state and a spin 1 state with the following 
result, 



i + ) = \l«m) + y-fPi), ;7:!! 

i.e., we are transforming the states l P\ and 3 Pi in the LS coupling scheme to the states 1 + ' 
and 1 + in the jj coupling scheme. 

In their discussions of the heavy quark fragmentation functions within the context of 



Heavy Quark Effective Theory, Jaffe and Randall JT9[ showed that a more natural variable 
to use is given by 

I _ f 

y = 1 , (75) 

r 

rather than the usual fragmentation variable z. Furthermore, these authors also showed that 
the heavy quark fragmentation function can be expanded as a power series in r, 

D(y) = ±a(y) + b(y)+0(r), (76) 

where a(y),b(y), etc. are functions of the variable y. The leading term a(y) is constrained 
by the heavy quark spin-flavor symmetry while all the higher order terms contain spin- 
flavor symmetry breaking effects. One can recast our results for the P-wave fragmentation 
functions derived in Section 3 in the above form. By carefully expanding the powers of r 
and (1 — fz), the leading contributions to the fragmentation functions are given by 

fl£W*) - ^^^(^-^ + 40^ + 96), (77) 

18 



- «^^T#(l/ 4 -% 3 + 8, 2 + 32), ,7* 



^- bc( 3 Pl) (y) - «^^-^(3 2/ 4 -V + 16 Z / 2 + 48), (79) 



243r y { 


5 


32 (y- 


i) 2 


243r 


! 


320 (y - 


I) 2 



D [ 1(V) - - ^^"''W + ^ + S). (81) 



Therefore, in the heavy quark limit, we have 



D 



^^i(y 4 -4y 3 + 8y 2 + 32), (82) 
Sir 



- «^^V + V + 8). (83) 
81r 

These imply the following spin counting ratios 



D- b l,i C (zp ){y) _ l 

z^- , M ~" 3 
• iy &_>& c ri+ / ) W 



6->6c(l+') 
/ 6_+6c(l+) 



^sS&W*/) 3 



J4) 



85) 



as expected from the heavy quark symmetry. 

Similarly, as already shown in Ref.||, the color-octet coefficients D^\^ c n S Az) and 

Di_^i c ns s dz) given by Eqn. (|62|) and (|63|) respectively also satisfy the heavy quark spin sym- 
metry in the limit of r — > 0. Thus, heavy quark spin symmetry is a powerful tool to check 
our tedious results obtained in previous Section. We note that the mixing effects and the 
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color-octet contributions were not considered in Ref. 0. Moreover, the results given there 
do not have the proper heavy quark limit given by Eqn.(fT^). 

Since our P-wave fragmentation functions derived in previous Section are consistent 
with the general QCD analysis of heavy quark fragmentation functions by Jaffe and Randall 
[ |19|j , they may be useful as phenomenological fragmentation functions for c — > D** and 
b — > B**, with the mass ratio r and the two overall factors o? s H\jm and a 2 s H' % /m treated as 
free parameters. 



5 Braaten-Levin spin counting rule 

The fragmentation function Di^u{z) for a parton i splitting into a hadron H is related to 
the distribution function /;/#(£) of finding the parton i inside the hadron H by analytic 



continuation 20, 21 



f l/H (x) = xD^ H (-) . (86) 

x 

The results of our perturbative S'-wave and P-wave fragmentation functions obtained in 
Ref.|| ^2[ and in this paper allow us to study the perturbative tail of the distribution 



functions of the heavy quark inside the heavy S-wave and P-wave mesons as well. From our 
explicit calculations, we see that fi/n{x) has a pole located at x — f. The pole is cut off 

by nonperturbative effects related to the formation of bound states of the (be) pair. This 
pole is of order 6 and 8 for the S'-wave and P-wave states respectively. In the general case 
of L-waves (L = 0, 1, 2, • • • or S-, P-, P-waves, • • •), we expect this pole is of order 6 + 2L. 
Therefore, we can expand f(x) as a Laurent series, 

f(x) = } + less singular terms , (87) 
(x — r) n 

with n = 6 + 2L for the general L-waves. 

An interesting result recently obtained by Braaten and Levin |23| states that the 
r-dependent coefficients a n (r) satisfy the simple spin counting rules. To be more specific, we 
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have 



for the S- waves, and 



aeCSo) : a 6 ( 3 5i) =1:3 



^Pi) : a 8 ( 3 P ) : a 8 ( 3 P 1 ) : a 8 ( 3 P 2 ) = 3:1:3:5 



(89) 



for the P-waves. 

From Eqns. (|62|) and ( p3] ) or using the results in Ref. |5|], we can calculate the two 
coefficients a 6 (}S ) and a 6 ( 3 Si) for the S-w&ve states to leading order in a s (up to an overall 
common factor of a 2 s H' & /m), 



aeCSo) 



a 6 ( 3 ^i) 



lim 

x— >r 



(x — r) xD\ 



So * 



X 



[rr] 



3 lim 

x— >r 



{x-f) 6 xDi%(- 



4 / 



(90) 
(91) 



Similarly, from Eqns.(]2"^), (|3"6]), Q3"7[), and (P8|), we obtain the coefficients a 8 ( 1 Pi) and a 8 ( 3 Pj) 
for the P-wave states to leading order in a s (up to an overall common factor of ot 2 s Hi/m), 



Cpi) 



as( 3 Pi] 



( 3 P, 



lim 

x—>r 



lim 

x—>r 



lim 

x— >r 



lim 



(a - ffxD^l (-) 
1 x 



(x-ffxDil(-) 

X 



(x-ffxDi 1 ^-) 



512 7 

rr 

81 V ; 



512 
243 1 



rr 



^7 



512 7 

rr) , 

81 K ' 



{x - ffxD^ 



x 



2560 
243 



(rr) 7 . 



(92) 
(93) 
(94) 
(95) 



Moreover, from Eqn. flBTD , we have 



lim 

x—>r 



[x - f) 8 xD { ± 



X 



(96) 



Therefore, our explicit results of the S'-wave and P-wave fragmentation functions do satisfy 
the Braaten-Levin spin counting rule. Notice that unlike the spin counting in the heavy quark 
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limit studied in the last Section which requires r — > 0, the Braaten-Levin spin counting rule 
holds for arbitrary values of r. 



6 Fragmentation functions for c — > h c and c — > XcJ {J = 

0,1,2) 

Our results derived in Sections 2 and 3 can be easily applied to the cases of P-wave charmo- 
nium and bottomonium by taking the equal mass limit of m c = rrib, i.e., r = f = 1/2. Simple 
formulas can be obtained in this limit and are given below for convenience. To be specific, 
we will consider the P-wave charmonium case. To leading order in v 2 , the fragmentation 
functions for a charm (or anti-charm) quark to fragment into various P-wave charmonium 
states consist of two terms: 

iW^o) = ^D^ 1Pi) (z,A)+3^^ D fl {z)t (97) 

D~**(z,IM>) = ^^ CC -(3P J) (^A) + (2J+1)%^^%3 51) (.), (98) 

where A is a factorization scale within the range m c v <C A < m c . Pi( c5 ) is related to the 
derivative of the radial wavefunction of the P-wave charmonium, P' P (0), according to JK]. [L3| 



Him « f M . (99) 



H's(cc) satisfies the renormalization group equation [10| , [14 



A ^X^( CC -)(A) = ^a s (A)H 1(cg) , (100) 



with the solution 

*WA) = ^8 (cc -)(Ao) + ^ lo s(^y) H litsSi , (101) 
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where (3q = 25/6 is the first coefficient in the beta function for QCD with 4 light flavors. If 
the factorization scale A is chosen to be much less than the charm quark mass m c , the above 
equation can be used to sum up large logarithms of m c /A. 

To leading order in a s , we have 



^ ) ci-cc( 1 Pi)( 2; ) 



V\2 



1 6aj(2m c ) z(l - z) 
81 (Z-zf 



x 



64 - 1282 + 176/ - 160/ + 140/ - 56/ + 9z t 



(102) 



16a 2 (2m c ) *(1 - 2) 
729 (2 - /) f 



x 



192 + 384^ + 528/ - 1376/ + 1060/ - 376/ + 59/ 



, (103) 



(i) 

C^CC^Pi) 



(*) = 



64a 2 (2m c ) z(l - zf 
243 (2 - zf 



x 



96 - 2882 + 496/ - 408/ + 202/ - 54/ + 7/ 



(104) 



D Hccpp 2 )( z ) 



128a*(2m c ) z(l - zf 
729 (2 - zf 



x 



48 - 192,2 + 480/ - 668/ + 541/ - 184/ + 23/ 



(105) 



The corresponding fragmentation probabilities = (Hi/m c ) Jq dzD^(z) are simply given 
by 



c— *h c 



p(l) 



p(l) 

c-*Xd 



p(l) 



8a 2 (2m e ) gi (ce) 
81 m c 

8a 2 (2m c ) jj"i (c c) 
729 m c 

64a 2 (2m c ) H 1(cE) 
243 m c 

32a 2 (2m c ) #i (c5) 
729 m r 



18107 
35 " 

119617 
35 

1151 



7 

54743 
35 



- 746 log 2 
-4926 log 2 

- 237 log 2 

- 2256 log 2 



(106) 
(107) 
(108) 
(109) 
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The color-octet coefficients can be obtained by making appropriate changes in the 
radial wavefunction and color factor to the results of the S"-wave fragmentation functions 
D c -*r) c ( z ) an d D c ^j/^{z) given in Ref . [ T2 ] , or simply by setting r = f = l/2in Eqns.(^) and 
(|63D. The results are 



D (8) (z) 



-Dcicc(3Si)( Z ) 



a 2 (2m c ) z(l - zf 
162 (2 - zf 

a 2 s (2m c ) z(l - zf 
162 (2 - zf 



48 + 8z 2 - 8z s + 3z 4 



16 - 32z + 72z 2 - 32z 6 + 5z 4 



(110) 



'III) 



The fragmentation probabilities P® = (H' 8 /m c ) J^dzD^(z) of the color-octet piece are 
given by 



p(8) 

^a-tcc^So) 



p(8) 

r c->cc( 3 Si) 



a 2 s (2m c ) H' 8[C5) (A) 
54 m c 

a](2m e ) K (cC) (A) 



^-371o g2 



162 



1189 



-57 log 2 



'1121 



(113) 



To avoid large logarithms of m c /A appearing in the color-singlet contributions, one should 
set A ~ m c in the matrix element HL^JA). 

Adding both color-singlet and color-octet contributions, the total fragmentation prob- 
abilities for c or c to split into the 4 P-wave charmonium state are given by 



c—fh c 



C^XcJ 



p(l) , qp(8) 

^c-tcc^Pi) "I" °- r c- + cE( 1 5 ) ' 



P 



(1) 



c->cc( 3 Pj) 



(114) 
(115) 



We note that in the equal mass case, mixings between the l P\ state and the 3 Pi state are 
not allowed by charge conjugation. Once produced by fragmentation, the P-wave Xcj states 
radiatively decay into J/ip and therefore contribute to the inclusive J/ip production cross 
section. 

Extensions of the formulas given in this Section to the P-wave states of bottomonium 
system are straightforward. 
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7 Discussions 

According to the results from potential model calculations || , the radial wavefunctions at the 
origin for the first two sets of P-wave (be) bound states are given by |P'2p(0)| 2 = 0.201 GeV 5 
and |i?3p(0)| 2 = 0.264 GeV 5 . The corresponding F 1(5c) (2) and H 1(hc) {3) are about 10 MeV 
and 14 MeV respectively, where we have used m c = 1.5 GeV and = 4.9 GeV. We will 
choose the factorization scale A = m. In the limit of m ^> Ao, the contribution to H' s ^ from 

the perturbative evolution dominates, and one can estimate HL^Jm) by setting a s (A ) ~ 
1 and neglecting the constant H' & ^JAq) in Eqn.(|). Taking a s (m) = 0.38, we estimate 
-^8(5e)( m ) ~ ^ ^ MeV and 1.8 MeV for the 2P and 3P states respectively. Numerically, 
the color-octet contributions are small compared with the color-singlet contributions. The 
mixing angle for the l P\ and 3 Pi states are also obtained in || with the results (cos 9, sin 9) 
= (0.999,0.030) and (0.957,0.290) for the 2P and 3P states respectively. The mixing angle 
is surprisingly small which implies the states constructed in the LS coupling scheme are 
actually very close to the physical eigenstates. The total fragmentation probabilities for b to 
split into the four 2P states are estimated to be 



Similarly, for the 3P states, we 



Pb-*bc(23p ) ~ 2.3 x 10 5 , 


(116) 


A^ C (21+')^ 4 - 4xl0 ~ 5 > 


(117) 


A-.M2i+)~4.8xlO- 5 , 


(118) 


Pb-*bc(23p 2 ) ~ 5.6 x 10~ 5 . 


(119) 


obtain 




FUic&*>) ~ 3 '° x ir5 ' 


(120) 


^c(31 + ')^ 8 - lxl0_5 ' 


(121) 


n^5c(3i + )~4.0xl0- 5 , 


(122) 


lUbcGPP*) « 7.4 x 10- 5 . 


(123) 



Since all the 2P states will decay 100% to the IS* pseudoscalar ground state B c , one 
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should add up all the probabilities from Eqns . ( |TT6| ) - ( |1 1 9| ) to give the total fragmentation 



probability for the B c production rate from the cascades of the 2P states, which is about 
1.7 x 10~ 4 . This is comparable to the probability 3.8 x 10~ 4 of the direct fragmentation of 
b — > B c , and is about 10 % of the lower bound 1.5 x 1CT 3 for the probability of b — > B c 
including the IS" and 2S states obtained in Ref. |J. Similarly, the total probability for the 
production of B c from the cascades of the 3P states is about 2.3 x 10~ 4 . Therefore the 2P 
and 3P states together account for a significant fraction of the inclusive production rate of 
the B c . While the 2P states are expected to lie below the BD flavor threshold, the 3P states 
may or may not lie below this threshold. If any one of the four 3P states lies above this flavor 
threshold, it will quickly dissociate into a pair of B and D mesons, and will not contribute 
to the inclusive B c production rate. The pr distributions and the total inclusive production 
cross sections for these P-wave states with realistic rapidity cut at the pp colliders are now 
under investigation |24| . The four 3-D states are also expected to lie below the BD threshold 



[[|. The corresponding P-wave fragmentation functions and probabilities have not yet been 
calculated. 

For the P-wave charmonium case, Pi( c5 ) has been phenomenologically determined 
to be ~ 15 MeV from the light hadronic decay rates of the xa and Xc2 [0- We set the 
factorization scale A = m c . Following Ref.|T2|, we will take H' 8 ^(m c ) pa 3 MeV. Numerically, 

the color-octet contributions are minuscule. The total fragmentation probabilities for c — > h c 
and c — > Xcj are estimated to be 



Pc- 


pa 1.8 x 1(T 5 , 


(124) 


Pc- 


>Xc0 Pa 2.4 x 10- 5 , 


(125) 


Pc- 


>Xcl Pa 2.8 x 10- 5 , 


(126) 


Pc- 


. Xc2 Pa 1.1 x 10- 5 . 


(127) 



The production of Xcj states contribute to the inclusive production rate of J / ip through the 
radiative decay Xcj —> J/ip + 7- Multiplying the fragmentation probabilities given above 
by the appropriate branching ratios of 0.7%, 27%, and 14%, we find that the probability 
of a J/ip inside a primary charm quark jet coming from the cascades of the Xcj states is 
approximately 0.9 x 10~ 5 . This is about an order of magnitude smaller than the probability 
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1.2 x 10 4 for the direct fragmentation of c — > J /if} obtained in Ref.p2| but still larger than 



the probability 3 x 10 for the direct fragmentation of g — > J /if) obtained in Ref. ||15||. The 
probability of finding a J /if} inside a primary gluon jet coming from the cascades of the Xcj 
states has also been estimated to be about 8 x 10 -5 ||12|| . The inclusive J /if) production rate 
at large transverse momentum region in pp collisions is expected to be dominated by the 
gluon fragmentation into Xcj followed by their radiative decays into J /if) [Q. Fragmentation 
contributions to the inclusive production cross sections of J/if), if)', and Xcj at pp colliders 
are now under investigation [EH . 

Fragmentation is the dominant mechanism for production of (be) mesons at large 
transverse momentum in high energy colliders. In this paper, we have calculated the frag- 
mentation functions for production of the P-wave (bc) mesons to leading order in a s . These 
P-wave fragmentation functions can be used to calculate the inclusive P-wave (be) mesons 
production cross sections at large transverse momentum. The production of the P-wave 
states account for about 20 % of the inclusive production rate of the ground state B c . The 
P-wave contributions to the inclusive production rate of B c should be smaller than those of 
the P-waves. Thus these P-wave fragmentation functions, together with the S'-wave fragmen- 
tation functions obtained in Ref. [Q, should allow an accurate calculation of the production 
rate of the B c meson at large transverse momentum in high energy colliders |24| . 
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Figure Caption 

1. Feynman diagram for b* — > bcc which contributes to the fragmentation of b into P-wave 
(be) bound states. The outgoing momenta are (1 — r)p + k, rp — k, and p' for the b, c, 
and c, respectively, k is the relative momentum of the b and c. 
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This figure "figl-l.png" is available in "png" format from: 
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